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ABSTRACT 

We investigate the influence of the star formation rate on the Tully-Fisher relation. We find 
that a simple model which combines the empirically-determined star-formation rate with the 
expected properties of galaxy halos provides a remarkably good fit to the absolute magnitude- 
rotation speed correlation. We find that the power-law nature, its slope, normalisation and 
scatter, are all readily accounted for if the Universe has a low density parameter, with or 
without a cosmological constant and disks are assembled at z ~ 1 — 1.5. Moreover, this 
agreement is found simultaneously in 4 wavebands. An Einstein-de Sitter Universe produces 
disks which are too faint unless the disks are assembled at z ~ 0.5. The scatter in the relation 
is due to a combination of the expected range of spin parameters of the halos and the range of 
formation redshifts. The source of the scatter opens up possibilities of a better galaxy distance 
indicator, if spectroscopic observations of globular clusters can be used to determine the halo 
rotation. 



1 INTRODUCTION 



The correlation between the luminosity and the rotation speed of 
spiral galaxies (L oc v a , with a ~ 3) has been used as a distance 
indicator for nearby galaxies (see e.g. Strauss & Willick 1995) 
since its discovery twenty years ago (Tully & Fisher 1977). In the 
very roughest terms, the correlation is explicable as a consequence 
of the virial theorem, but the fact that the slope of the relation is 
dependent on the waveband used suggests that the picture is at best 
slightly more complicated. There are clearly at least two potentially 
relevant ingredients in a model to explain the Tully-Fisher relation, 
these being the properties of the mass distribution (the halo) and the 
stellar properties of the baryonic disk material. Mo, Mao & White 
(1998) assume the disk mass-to-light ratio is constant, and therefore 
argue that it is the halo properties which control the Tully-Fisher re- 
lation. On the other hand Silk (1997) argues that self-regulated star 
formation is the controlling influence, and the Tully-Fisher relation 
has nothing to do with the halo. The picture may be complicated 
still further by effects such as energy feedback from stars, and fresh 
infall of gas over time, motivating an approach where all physical 
processes are modelled with complex computer codes (e.g. Stein- 
metz & Navarro 1998). In this paper, we take a simpler comple- 
mentary approach, which combines halo properties with empirical 
star formation properties. This approach has the advantage that it 
offers a possibility of improving the correlation by measuring fur- 
ther physical parameters; this is the main motivation for taking this 
analytic approach. It also appears to account for the relation ex- 
tremely well. We model the halo with the theoretically known spin 
properties of collapsed objects, using the isothermal sphere for the 
density profile, partly for simplicity, and partly because the issue of 
the dark matter profile in collapsed haloes is not yet settled (Moore 
et al. 1998). 

We model the star formation entirely empirically, using the 
Schmidt law relating star formation rate to disk surface density 
(Kennicutt 1998). With this simple model, we are able to obtain an- 
alytic expressions for the star formation rate as a function of time. 



The spectrophotometric evolution cannot be done accurately in an 
analytic way, so for this step we use a synthetic stellar population 
numerical code (Jimenez et al. 1998; Jimenez et al. 1999). Rather 
surprisingly, this simple model is able to account extremely well 
for the Tully-Fisher relation. We are able to reproduce the slope, 
normalisation and scatter of the relation, simultaneously in 4 wave- 
bands, with very little freedom. The idea of the spin parameter in- 
fluencing the galaxy properties is of course not new; it has been 
standard since Fall & Efstathiou (1980), and applications have been 
made for disk galaxies for example by Kashlinsky (1982), Dalcan- 
ton, Spergel & Summers (1997), Jimenez et al. (1997), Jimenez 
et al. (1998), as well as Mo, Mao & White (1998); van den Bosch 
(1998). The new element in the analysis here is the combination 
with the empirical dependence of star formation rate on surface 
density. 



2 HALO AND STAR-FORMATION MODEL 

We follow the notation and method of Mo, Mao & White (1998) 
in the modelling of the halo by an isothermal sphere, characterised 
by its mass M, and circular velocity V c - These are related via the 
Hubble constant at the redshift z of formation H(z) 



M = 



V? 



WGH(z) 



(1) 



Note that V c is not the actual circular velocity of the halo - it is the 
circular velocity required for centrifugal support in the potential of 
the halo. In terms of the present Hubble constant and the density 
parame- 
ter in non-relativistic matter and cosmological constant, H(z) = 

H [Hao + (1 - Qao - fio)(l + zf + Q (l + z) 3 ] 1/2 . 
The disk mass we assume to be Md = ra^M, where the frac- 
tion m,d = fib /Oq and the baryon density is set by nucleosynthesis 
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(Walker et al. 1991), Q. b = 0.015ft' 2 , where ft = H /WO km s" 1 
Mpc -1 is taken to be 0.65 throughout. 

We assume the disk settles with an exponential profile 
E(7?) = Eo exp(— R/Rd), with central surface density Eo and 
scale length Rd, in which case Md = 27rEo_R^. The disk scale 
length is related to the spin parameter A of the halo by (Mo, Mao 
& White 1998) 



Rd 



WV2H(z) 



(2) 



where we have assumed that the disk and halo have the same spe- 
cific angular momentum (Mestel 1963). A = J|S| 1/2 G _1 M" 5/2 , 
where E is the total energy of the halo, and J its angular momen- 
tum. This sets the initial gas surface density, which we assume is 
accreted on a timescale which is short compared with a Hubble 
time. We also ignore any gas returned to the ISM by stars, or late 
infall of fresh gas. 

We assume that the star formation rate is set by the empiri- 
cal Schmidt law (Kennicutt 1998), dependent only on the local gas 
surface density E 9 , * S fr = BEg 1 ' 4 * 015 : 



*sfr = (2.5±0.7)xl0 



-(- 



\ M Q pc~ 



M @ y- 1 k P c-' 2 .(3) 



The time-evolution of the gas surface density is therefore given in 
terms of the initial gas surface density E 9 o by 

-2.5 



J 90 



+ 0.4Bt) 



(4) 



where B — 9.5 x 10 in SI units and t is measured from the 
assembly of the disk. From this and the star formation law, we can 
integrate over the disk to compute the total star formation rate: 

M*W = ™y° 4 f .////,-xp.V • a) ' (5) 



(in kg s 1 ), where a(t, Eo) = 0.4f?tEo' 4 . This may be written in 
terms of a generalised hypergeometric function (e.g. Gradshteyn & 
Ryzhik 1980): 



M,(t) 



507rBfl 2 E^ 4 
49 



jF 2 (3.5, 3.5, 3.5; 4.5, 4.5; -a) 



(6) 



A similar integral gives the remaining gas mass as a function of 
time as 



M g (t) = 27ri?iE 3F2 (2.5, 2.5, 2.5; 3.5, 3.5; -a) . 



(7) 



In more astronomical units: 

M,(t) = 247ft 



o.e (J2!±\ ( 

V0.05/ \l 



250kms~ 1 
3 F 2 (3.5, 3.5, 3.5; 4.5, 4.5; -a) M^y' 1 



(0.O5) 



x 
(8) 



and 



1.06ft" 



Gyr) ( 



(<: ) V 250A-///.S 



h z = H(z) /WO km s _1 Mpc -1 . We therefore have a complete 
description for the star formation rate as a function of time, for any 
cosmology and any halo parameters V c , z and A. The relative star 
formation rate and remaining gas fraction are shown in Fig. 1. The 
gas fractions remaining depend sensitively on cosmology and cir- 
cular velocity but are typically in the range 1 — 10% (c.f. Roberts 
& Haynes (1994)). The star formation rate is then fed into a sophis- 
ticated spectrophotometric stellar evolution code, which has been 
described elsewhere (Jimenez et al. 1998; Jimenez et al. 1999). 



remaining gos 
star formation rate 




VO.05/ VO.05/ v ' 



Figure 1. The star formation rate and remaining gas fraction, normalised to 
unity at the disk formation time, as a function of the parameter a <x t. 



Broad-band magnitudes inB,R,I and K' are obtained from the re- 
sulting spectra. Note we have chosen a particular synthetic popula- 
tion code; other codes, such as those discussed in Chariot, Worthey 
& Bressan (1996), agree to within 0.2 magnitudes. This leads to a 
small shift in normalisation, but does not add dispersion. We as- 
sume a solar metallicity throughout. Including chemical evolution 
would lead to higher spin systems being bluer, and lower spin sys- 
tems being redder than we have computed, but the effect is confined 
toB. 



3 RESULTS 

We show the results as graphs of absolute magnitude against circu- 
lar velocity, in B, R, I and K', for disk formation redshifts z=l,2,3 
and spin parameters A=0.025, 0.05, 0.1. These spin parameters rep- 
resent the 10, 50 and 90 percentile points for a gaussian field for any 
reasonable power spectrum (e.g. Warren et al. 1992). Disks are sub- 
ject to disruption by major merging events, and there are arguments 
that present-day disks may have been assembled at z ~ 1 — 2 (Weil, 
Eke & Efstathiou 1998), and there are arguments from disk sizes 
that argue for even lower assembly redshifts (Mo, Mao & White 
1998). In Fig. 2, we show the results for a flat cosmological model 
with h — 0.65, Qo = 0.3, fUo = 0.7. Superimposed is the data 
from the careful study of spiral galaxies in the loose clusters of 
Ursa Major and Pisces (Tully et al. 1998). The inclination-corrected 
FWHM of the lines has been converted there to W % R , which approx- 
imates to twice the circular velocity. From Figure 2 and also for the 
model of Figure 4, we find the following: 

• The T-F relation is a good power-law, with about the right 
slope 

• The T-F normalisation is reproduced for plausible formation 
redshifts 

• The scatter is comparable to observation, if disk formation 
redshifts are not too widely distributed 

Moreover, these successes apply simultaneously to all four colours 
(especially for R, I and K' , which sample the old population). The 
blue data are a little more ragged; nevertheless our model provides 
a good fit leaving little room for any enhanced star formation due 
to late infall of gas. It is worth noting that the predicted level of star 
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Figure 2. The absolute magnitude (in 4 bands) vs. circular velocity for disk 
galaxies in a flat, A-dominated universe, with f2o = 0.3, f2y\o = 0.7, and 
Ho = 65 km s" 1 Mpc" 1 . The solid, dotted and dashed lines correspond 
to disk formation redshifts of 1,1.5 and 2 respectively. For each formation 
redshift, the 3 lines show halos with spin parameters A = 0.025, 0.05 and 
0.1 (bottom-to-top). These spin parameters are the 10, 50 and 90 percentile 
points for the spin distribution from gaussian initial conditions. The data are 
drawn from Ursa Major and Pisces clusters (Tully et al. 1998). Also plotted 
are the near-dark galaxies NGC2915 (diamond) and DD0154 (triangle). 



formation agrees well with that found by Madau (1997); Steidel 
et al. (1998) since our model has typical initial SFR of 10-100 M Q 
y _1 with a highest possible value of 1000 M© y _1 for the systems 
with the highest V c and lowest A. 

The magnitude is not necessarily a monotonic function of spin 
parameter, at fixed V c . Low-spin systems forming at given z will 
have high surface densities and high initial star formation rates. At 
first they will be brighter than the high-spin systems. However, the 
timescale over which they use up most of their gas is short, and the 
fading of stars over time may make them fainter than the high-spin 
systems at late times. 

Fig. 3 shows results for an Einstein-de Sitter cosmology. This 
model forms disks with very high surface densities, which have 
high initial star formation rates, but which fade excessively unless 
their formation redshift is low (~ 0.5). One could argue that we 
have a small amount of freedom in moving the curves left and right, 
since the observed relation uses the FWHM W of the lines, and we 
have assumed W = 214. This is probably most appropriate for 
the isothermal sphere model we assume, but one might allow some 



flexibility here. We have also investigated an open universe with 
zero cosmological constant; it fits almost as well as the A model. 

In Fig. 4, we show the flat model, but we incorporate the peak 
height-spin parameter anticorrelation (approximately A oc v~ 1A ) 
claimed from analytical studies by Heavens & Peacock (1988) (see 
also Steinmetz & Bartelmann 1994). There is some support for 
this from the numerical simulations of Ueda et al. (1994). Lemson 
& Kauffmann (1999) see no anticorrelation with density, but they 
present results only for density in a rather large sphere (10 
Mpc). Since the density on galaxy scales (filter length ~ 0.5h~ 
Mpc) is only weakly correlated with density on such large scales, 
it is hard to see how any small-scale correlation could be appar- 
ent in such a study. Assuming that the small-scale anticorrelation 
is present makes little difference to the Tully-Fisher relation, for a 
reasonable mean formation redshift, which we take to be 1.5 (Note 
that this refers to the mean collapse redshift of the haloes; the disks 
observable today may be younger and have a lower assembly red- 
shift). Early-forming high peaks have systematically lower spin, 
which leads to higher surface density, brighter disks. The predicted 
scatter in the relation depends on the probability distribution of the 
disk assembly redshift. We can estimate a lower limit to the scat- 
ter by assuming a single formation redshift. For 2 = 1— 1.5, the 
scatter is approximately 0.3 magnitudes in K and 0.7 in B; the ob- 
served values in Ursa Major and Pisces vary from 0.42 in R, I and 
K' to 0.55 in B, although a sample of 12 clusters shows a smaller 
scatter of 0.35 in / (Tully 1998). Inspection of Fig. 4 shows that 
a range of formation redshifts can give good agreement with the 
observed scatter. Note that the uncertainty in the Kennicutt fit to 
the Schmidt law (equation [3]) makes relatively little impact on the 
results; changing the normalisation of the star formation rate by 30 
per cent either way alters the magnitudes by about 0.15. Further we 
note that the anticorrelation would ease the problem of disks being 
too small compared with observation (Mo, Mao & White 1998), 
since the later-forming haloes have typically larger spins than aver- 
age. Late-forming haloes forming at redshift 1 (when the mean is 
1.5) are 40% larger if the anticorrelation holds. 

The thick lines in Fig. 4 show the theoretical limit for galaxies 
for forming stars according to the Toomre (1964) stability criterion 
(see also Kennicutt 1989). As discussed in Jimenez et al. (1997) 
halos with V c < 40 km s~ fail to form stars for any value of A, 
while for halos with V c > 120 km s _1 stars are formed in the 
disk for any value of A as shown in Fig. 4 of Jimenez et al. (1997). 
The absence of faint discs is also present in other samples (see e.g. 
Strauss & Willick 1995), although observational selection effects 
may be the explanation for current data (Tully, priv. comm.). To 
illustrate this we have also plotted the near-dark galaxies NGC2915 
(diamond) and DD0154 (triangle) (see Meurer et al. 1996). These 
systems fall quite well in the scenario described in Jimenez et al. 
(1997) since only the central parts of these galaxies have formed 
stars, and have a dark HI disk. The reason why they fall below the 
theoretical TF law is due to the fact that most of its HI mass is not 
in stars (they have mass-to-light ratios of about 80), thus they are 
under-luminous. 

Fig. 5 shows the predicted Tully-Fisher relation for disks ob- 
served at redshift 1. If the disk assembly is also at redshift 1, then 
the disks are extremely bright in all wavebands. For disk assembly 
at redshifts 1.5 and 2, the Tully-Fisher relation is generally a little 
brighter, the degree of brightening depending in detail on circular 
velocity, band and formation redshift. For massive systems forming 
at redshift 1.5, the systems are about 0.5 — 1.5 magnitude brighter 
in blue, somewhat larger than the observations of Vogt et al. 1997, 
although the samples are not large enough as yet for the situation 
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Figure 3. As Fig. 2, but for an Einstein-de Sitter universe and formation 
redshifts of 0.5 (solid line), 1.0 (dotted line) and 1.5 (dashed line). 



to be clear (Rix et al. 1997, Hudson et al. 1998, Simard & Pritchet 
1998). 

The model predicts some dependence of the Tully-Fisher rela- 
tion on the surface brightness of the galaxy, in the sense that higher 
surface brightness disks should have higher luminosities for given 
rotation speed. The effect is not large, as witnessed by the tightness 
of the predicted relations in Fig. 4, and there is also some overlap in 
the surface brighnesses at fixed M, V c arising from a spread in disk 
assembly redshifts. Thus the coincidence of TF relations of high 
and low surface brightness disks (Sprayberry et al. 1995; Zwaan 
et al. 1995) presents no great difficulty for the model. More serious 
is the bimodal surface brightness distribution of disks observed by 
Tully & Verheijen (1997) (but see Mihos, Spaans & McGaugh 1998 
for an alternative view), which does not appear naturally in this 
model. The authors suggest that the high surface brightness disks 
may be self-gravitating, whereas the LSBs are dark-matter domi- 
nated. We cannot exclude such a possibility, as our model assumes 
a given halo potential. If the passage to self-gravitation means only 
more early star formation (from higher surface densities), then the 
results we present are virtually unchanged: at these wavebands, the 
distinction between a burst of star formation and the ongoing trickle 
which our model predicts is negligible for the high surface density 
disks. Observations at U would be required to distinguish the his- 
tories. 
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Figure 4. As Fig. 2, but including the halo peak height-spin parameter an- 
ticorrelation predicted by Heavens & Peacock (1988). The thick solid lines 
shows the expected theoretical threshold for halos that fail to form stars 
(Jimenez et al. 1997), for formation redshift 1.5. System below the line 
will fail to form any significant population of stars and thus will be 'dark'. 
Comfortably enough, no galaxies are observed below the threshold in any 4 
bands, except for system with stars confined to the central regions with few 
stars in the disk, such as NGC2915 (diamond) and DD0154 (triangle). Also 
plotted as open squares are the data for LSBs from Zwaan et al. (1995), this 
illustrates the point that LSB and HSB galaxies lie in the same TF law 



4 SUMMARY AND PROSPECTS FOR A BETTER 
DISTANCE INDICATOR 

We have combined the theoretical properties of dark matter halos 
with the empirically-determined dependence of star formation rate 
on disk surface density to predict a correlation between luminos- 
ity and circular velocity which is compared with observation. Re- 
markably good agreement is found, in 4 wavebands, with the re- 
cent study of Tully et al. (1998), for a flat universe dominated by 
a cosmological constant if the disk assembly redshift is 1 — 1.5, 
in agreement with conclusions from Weil, Eke & Efstathiou (1998) 
which are based independently on disk survival rates. An Einstein- 
de Sitter universe can be reconciled if the disk formation redshift 
is low: z ~ 0.5. This explanation of the Tully-Fisher relation falls 
between the halo-determined extreme of Mo, Mao & White (1998) 
and the halo-independent model of Silk (1997), however, it is only 
a partial explanation, as it relies on the empirical Schmidt law (see 
Silk 1997 for some relevant discussion). This study suggests that 
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Figure 5. Predicted Tully-Fisher relation at z = 1 for the A universe of 
Fig. 2 and formation redshifts of 1.5 (dotted line) and 2 (dashed line). 



there are two parameters which lead to scatter in the Tully-Fisher 
relation. These are the spin parameter of the halo, and the redshift 
at which the disk is accumulated. The age of the disc appears to be 
less important than the redshift, which controls the initial surface 
density of the disk and hence the star formation rate. The general 
effect of the formation redshift and the spin parameter is almost in- 
evitable, and is likely to be more robust than the specific model cal- 
culations presented here. The relative influence of formation red- 
shift and spin parameter is model-dependent. In particular, if, as has 
been claimed (Heavens & Peacock 1988; Ueda et al. 1994), there 
is an anticorrelation between peak overdensity and spin parameter, 
then the predicted Tully-Fisher relation is changed only slightly. 
The scatter might be narrowed by measurement of the age and spin 
parameter. The latter is a difficult quantity to measure, but could be 
attempted by looking at the systematic velocity of halo stars which 
are not formed dissipatively in the disk. Individual stars would be 
impossible to measure, especially with a bright disk present, but it 
is possible that for nearby galaxies the integrated light from globu- 
lar clusters could be used. 
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